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I. INTRODUCTION 

The study of topological insulators has become an im- 
portant research focus in condensed matter physics--'?. 
Historically, the research of topological band insulators 
started from the two dimensional (2D) quantum Hall ef- 
fect. Landau level (LL) quantization gives rise to nontriv- 
ial band topology characterized by integer-valued Chern 
numbers .-iii^ In fact, LLs are not the only possibility for 
realizing topological band structures. Quantum anoma- 
lous Hall band insulators with the regular Bloch-wave 
structure are in the same topological class as 2D LL sys- 
tems in magnetic fields^. Later developments generalize 
the anomalous Hall insulators to time-reversal (TR) in- 
variant systems in both two and three dimensions i^^— 
This is a new class of topological band insulators with TR 
symmetry which are characterized by the Z2 index. Ex- 
perimentally, the most obvious signatures of band topol- 
ogy appear on open boundaries, in which they exhibit he- 
lical edge or surface states. Various 2D and 3D materials 
are identified as topological insulators, and their stable 
helical boundary modes have been detectecl?lS-^22. Fur- 
thermore, systematic classifications have been performed 
in topological insulators and superconductors in all the 
spatial dimensions, which contain ten different universal 
classea^ii^. 

Although the current research is mostly interested in 
topological insulators with Bloch-wave band structures, 
the advantages of LLs make them appealing for further 
studies. We use the terminology of LLs here in the fol- 
lowing general sense not just for the usual 2D LLs in 
magnetic fields: topological single-particle level structures 
labeled by angular momentum quantum numbers with 
flat or nearly flat spectra. On open boundaries, LL sys- 
tems develop gapless surface or edge modes which are 
robust against disorders. For example, in the 2D quan- 
tum Hall LL systems, chiral edge states are responsible 
for quantized charge transport. For the 2D LL based 
quantum spin Hall systems, helical edge modes are ro- 



bust against TR-invariant disorders^. Similar topological 
properties are expected for even high-dimensional LL sys- 
tems, which exhibit stable gapless surface modes. For the 
usual 2D LLs, the symmetric gauge is used in which an- 
gular momentum is conserved. We do not use the Landau 
gauge because it does not maintain rotational symmetry 
explicitly. LL wavefunctions are simple and explicit, and 
their elegant analytical properties nicely provide a plat- 
form for further study of topological many-body states 
in high dimensions. 

Generalizing LLs to high dimensions started by Zhang 
and Hu^'^ on the compact sphere by coupling 
large spin fermions to the SU(2) magnetic monopole, 
where fermion spin scales with the radius as . 
Later on various generalizations to other manifold were 
developedi^"— Two of the authors have generalized the 
LLs of non-relativistic fermions to arbitrary dimensional 
flat space i?'^ —. The general strategy is very simple: the 
harmonic oscillator plus spin-orbit (SO) coupling LijTij, 
where Lij and are the orbital and spin angular mo- 
menta in a general dimension. Reducing back to two 
dimensions, it becomes the quantum spin Hall Hamilto- 
nian in which each spin component exhibits the usual 2D 
LLs in the symmetric gauge, but the chiralities are oppo- 
site for two spin components®. For a concrete example, 
say, in three dimensions, each LL contributes a branch of 
helical Dirac surface modes at the open boundary, thus 
its topology belong to the 2'2-class. Furthermore, LLs 
have also been constructed to arbitrary dimensional flat 
spaces for relativistic fermions'^'', which is a square root 
problem of the above non-relativistic cases. It is a gener- 
alization of the quantum Hall effect in graphene^"— to 
high dimensional systems with the full rotational sym- 
metry. This construction can also be viewed as a gener- 
alization of the Dirac equation from momentum space to 
phase space by replacing the momentum operator with 
the creation and annihilation operators of phonons. The 
zero-energy LL is a branch of half-fermion modes. When 
it is empty or fully occupied, fermions are pumped from 



2 



the vacuum, a generalization of parity anomaly^— to 
high dimensions. 

In this article, we study another class of isotropic LLs 
with TR symmetry but breaking parity in two and three 
dimensions, which can also be straightforwardly general- 
ized to arbitrary dimensions. The Hamiltonians are again 
harmonic oscillator plus SO couplings, but here the SO 
coupling is the coupling between spin and linear momen- 
tum, not orbital momentum. In 2D, it is simply the stan- 
dard Rashba SO coupling, and in 3D it is the a-p-type SO 
coupling. In both cases, parity is broken. The strong SO 
coupling provides the projection of the low energy Hilbert 
space composed of states with the proper helicity. The 
radial quantization from the harmonic potential further 
generates gaps between LLs. The SO coupling strongly 
suppresses the dispersion with respect to the angular mo- 
mentum within each LL. In two and three dimensions, 
they exhibit gapless helical boundary modes which are 
stable against TR-invariant perturbations, thus they be- 
long to the Z2 topological class. In fact, parent Hamil- 
tonians, whose first LL wavefunctions are obtained ana- 
lytically and whose spectra are exactly flat, can be con- 
structed by the dimensional reduction method from the 
hi^h-dimensional LL Hamiltonians constructed in Ref. 

This paper is organized as follows. The study of 
isotropic and TR-invariant LLs with parity breaking is 
presented in Sect. [Ill The generalization to three dimen- 
sions is given in Sec. IIIII The experimental realization of 
the 3D Rashba- like a ■ p -type SO coupling is performed 
in Sec. IIVI Conclusions and outlook are summarized in 
Sec. El 



II. TWO-DIMENSIONAL SPIN-ORBIT 
COUPLED LANDAU LEVELS WITH 
HARMONIC POTENTIAL 

In this section, we consider the Hamiltonian of Rashba 
SO coupling combined with a harmonic potential 

H2D = - + -mw^r^ - Xi-ihVxUy + ihVyax), (1) 

where u> is the trapping frequency; A is the SO coupling 
strength with the unit of velocity. Equation ([T]) is in- 
variant under the 5*0(2) rotation and the vertical-plane 
mirror reflection. In other words, the system enjoys C^oo 
symmetry. Equation ([I} also satisfies the TR symme- 
try of fermions, i.e., T = ia2K, with — —1 and K 
the complex conjugation. However, parity symmetry is 
broken explicitly by the Rashba term. 

Equation[T]can be realized in solid-state quantum wells 
and ultra-cold atomic traps. Rashba SO coupling due to 
inversion symmetry breaking at 2D interfaces has been 
studied extensively in the condensed matter literaturej^ 
its energy scale can reach very large valuesi^ Further- 
more, Wigner crystallization in the presence of Rashba 



SO coupling has been studied.^ In the context of ul- 
tra cold atoms, Bose-Einstein condensation with Rashba 
SO coupling plus harmonic potential was studied by one 
of the authors and Mondragon in Ref. in which 

the spontaneous generation of a half-quantum vortex 
is found. Later, there was great experimental progress 
in generating a synthetic gauge field from light-atom 
interactiouj^? which inspired a great deal of theoretical 
interest 4^"— 

A. Energy spectra 

In a homogeneous system with Rashba SO coupling, 
i.e., w = in Eq. ([T]), the single-particle states ip±{k) 
are eigenstates of the helicity operator a ■ {k x z) with 
eigenvalues ±1, respectively. The spectra for these two 
branches are e±(fc) = h^{k =F fco)^/(2m), and the low- 
est energy states are located around a ring with radius 
fco = mX/h in momentum space. Such a system has two 
length scales: the characteristic length of the harmonic 

trap It = \J and the SO length scale Igo = l/^o- 
The dimensionless parameter a = It 1 1 so describes the 
SO coupling strength with respect to the harmonic po- 
tential. 

As presented in Ref. [41)] for the case of strong SO 
coupling, i.e., 1, the physics picture is mostly clear 
in momentum representation. The lowest energy states 
are reorganized from the plane- wave states 'ip+{k) with 
k near the SO ring. Energetically, these states are sepa- 
rated from the opposite-helicity ones V'- [k) at the order 
of Ego = hk^X = a^Etp, where Etj, = Huj is the scale 
of the trapping energy. As shown below, the band gap 
in such a system is at the scale of Etp. Since a ^ 1, 
we can safely project out the negative helicity states 
ijj-^k). After the projection, the harmonic potential in 
momentum representation becomes Laplacian coupled to 
a Berry connection Ak as 

Vtp = ^uj\i^k~Ak)\ (2) 

which drives particle moving around the ring with a mo- 
ment of inertial / = Mkk'^; Mk — h^/[muj'^) is the ef- 
fective mass in momentum representation. The Berry 
connection is defined as 

Ak = i{'4)k+\^ k\'tpk+) = ^efc, (3) 

where |V'fe+) is the lower branch eigenstate with momen- 
tum k. It is well known that for the Rashba Hamiltonian, 
the Berry connection Ak gives rise to a tt flux at fc = (0, 0) 

but without Berry curvature at A: ^ 0.— This is because 
a two-component spinor after a 360° rotation docs not 
come back to itself but acquires a minus sign. 

The crucial effect of the tt flux in momentum space 
is that the angular momentum eigenvalues become half- 
integers as jz = m + The angular dispersion of the 
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spectra becomes Eagi{jz) = = {jl/2a'^)Etp. On 

the other hand, the radial potential in momentum repre- 
sentation is V{k) = ^Mkijj'^{k — fco)^ for positive-helicity 
states. For states with energies much lower than Eso, we 
approximate V{k) as harmonic potential, thus the radial 
quantization is Erad{nr) — {fir + h)Etp up to a constant. 
The same dispersion structure was also noted in recent 
works^"— , which show 

i^«...-(-. + ^-Y + |j)i^*P, (4) 

where the zero point energy is restored here. Since a ^ 
1, we treat rir as a band index and jz as a good quantum 
number for labeling states inside each band. 

B. Dimensional reduction from the 3D Landau 
level Hamiltonian 

Equation ([1]) not only can be introduced from the solid- 
state and cold atom physics contexts, but also can be 
viewed as a result of dimensional reduction from a 3D LL 
Hamiltonian [Eq. ([5])] proposed by the authors in Ref. 
i29i . This method builds up the connection of two topo- 
logical Hamiltonians in three dimensions with inversion 
symmetry and two dimensons with inversion symmetry 
breaking. The resultant 2D Hamiltonian Eq. (O exhibits 
the same physics that Eq. ^ does for eigenstates with 
jz < a in the case of a >• 1. The advantage of Eq. ([7]) is 
that its lowest LL wavefunctions are analytically solvable 
and their spectra are flat. 

Just like the usual 2D LL Hamiltonian in the symmet- 
ric gauge, which is equivalent to a 2D harmonic oscillator 
plus the orbital Zeeman term, the 3D LL Hamiltonian is 
as simple as a 3D harmonic potential plus SO coupling^^ 

H3D,LL = 7: h -muj'^r^ - ujL ■ a, (5) 

2m 2 

which possesses 3D rotational symmetry and TR sym- 
metry. Its eigen-solutions are classified into positive- and 
negative-helicity channels according to the eigenvalues of 
a ■ L = Ih or —{l + l)h, respectively. In the positive 
(negative)-helicity channel, the total angular momentum 
j± = (/ ± ^)h. The spectra in the positive-helicity chan- 
nel, En^^i = {2nr + are dispersionless with respect 
to the value of j+, thus these states are LLs. In the pres- 
ence of an open boundary, each filled LL contributes a 
branch of helical Dirac Fermi surface described as 

Hsf ^ Vf{(J X p) ■ ir - H, (6) 



where Br is the local normal direction of the surface, Vf 
the Fermi velocity, and fi the chemical potential. The 
stability of surface states under TR-invariant perturba- 
tions are characterized by the Z2 topological index. 

Now let us perform the dimension reduction on Eq. ([5]) 
by cutting a 2D off-centered plane perpendicular to the 
z-axis with the interception zq. Within this 2D plane of 
z = zo, Eq. ([5]) reduces to 

El2D,re = H2D-'^LzCrz- (7) 

The first term is just Eq. ([T|) with Rashba SO strength 
A = ujzq, and the 2D harmonic trap frequency is the same 
as the coefficient of the Lz(Tz term. The dimensionless 
parameter a — It/Iso — \zo\/It- If zq — 0, Rashba SO 
coupling vanishes. In this case, Eq. ([7]) becomes the 
2D quantum spin Hall Hamiltonian proposed in Ref. |^, 
which is a double copy of the usual 2D LL with opposite 
chiralities for spin-up and -down components. At zq ^ 0, 
Rashba coupling appears which breaks the conservation 
of (Jz- 

Two of the authors found the lowest LL solutions for 
Eq. ([SJ, whose center is shifted from the origin to Tc — 
(0, 0, 2;o) in Ref. 11^ These states do not keep j conserved 
but do maintain jz as a good quantum number as 

X { , tiY^ , V (8) 
\^ -sgn(zo)e'''Jm+i(fco/o) / 



where p = \/ + y^; jz = m + ^; fco = zq/^tJ and is 
the azimuthal angle around the z axis. The ip3D.j^,zo^s 
form a complete set of the lowest LL wave functions, but 
they are nonorthogonal if their j^'s are the same. By 
setting z — Zq in the above wavefunctions, we define the 
2D reduced wave functions as 

^^"-^^'^)=' ^ l-sgn(.o)e^("+^)^J„.+i(W) j- 

(9) 

Noticing that dzip3D,j^,zo\z=zo — 0, it is straightforward 
to check that i/'2D.j^ 's are solutions for the lowest LLs for 
the 2D reduced Hamiltonian in Eq. ([7]) as 

2 

H2D,re ^2D.j, = ~ ^ V'2D,i,- (10) 

The TR partner of Eq. © can be written as 



2D, 



sgn(zo)e-*("+i)'^J™+i(fcop) 



sgn(zoJe 



-i{m+l)(l) J 



(m+1) 



~sgn(zo)e-""^J_™(fcop) 



(11) 
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FIG. 1: (Color online) Energy dispersions of the solutions 
for the first four LLs to the 2D reduced Hamiltonian Eq. ((7} 
(solid lines), and those for (Eq. [1} (dashed lines). The value of 
a = It /ho ~ 35. The lowest LLs of Eq. Q are dispersionless 
with respect to j^. Please note that the overall shift of the 
zero-point energy difference |fia; is performed for the spectra 
of Eq. IH) for a better illustration. 

C. Relation between Eq. ([T]) and Eq. ([7] 

) 

The difference between the two Hamiltonians, Eq. ([7]) 
and (Eq. [1]), is the L^az term. Its effect depends on 
the distance p from the center. We are interested in the 
case of |zo| ^ It, «-e., a » 1. Let us first consider 
the lowest LL. With small values of jz, i.e., m < a, 
Jm{kop) and Jm+i(fcop) already decay before reaching 
the characteristic length Ix of the Gaussian factor. We 
approximate their classic orbital radii as the locations 
of the maxima of Bessel functions, which are roughly 
Pcj^ ~ ^It < It- In this regime, the effect of LzCTz 
compared to the Rashba part is a small perturbation, of 
the order of ojpcj^/X = PcjJzq <^ 1. Thus, these two 
Hamiltonians share the same physics. On the other hand, 
let us consider the case of very large values of jz, say, 
m ^ a^. The Bessel function behaves like p™ or 
at < p < ^It- The classic orbit radii are just pcj^ ~ 
y/mlT- The physics of Eq. (O in this regime is dominated 
by the Lz<Jz term and, thus, is the same as that of 2D 
quantum spin Hall LL wave functions. However, for Eq. 
([T]), the projection to the sub-Hilbert space spanned by 
Tp+ik) is not valid. Its eigenstates in this regime cannot 
be viewed as LLs anymore. For intermediate values of 
jz, i.e., a < m < , the physics is a crossover between 
the above two limits. 

For higher LLs of Eqs. ^ and ^ , we expect that their 
wave functions can be approximated by a form of Eq. Q 
by multiplying a polynomial of p at the n^-th power. As 
a result, the physics is similar to what is analyzed in 
the previous paragraph. At small values of jz < a, the 
energy gap is quantized in terms of the unit of Etp — fiw 
as in Eq. ^ for both Hamiltonians. At very large values 
of ^2 ^ a^, the LLs of Eq. ([7]) become flat again and the 
quantization gap is at 2Etp — 2hjj. 

We perform the numerical calculation of the energy 
levels of the reduced 2D Hamiltonian Eq. ([7]), as plotted 



in Fig. [T] The numerically calculated spectra of Eq. [H 
which were plotted in Refs. 113 andlH, are also presented 
for comparison. Only the spectra of jz > are plotted, 
and those of jz < are degenerate with their partners 
by the TR transformation which flips the sign of jz ■ The 
lowest LL of Eq. ([7]) is flat as expected, while higher LLs 
are weakly dispersive which is hardly observable for the 
range of jz presented. The LLs of Eq. (JT]) are dispersive 
with the dependence on jz shown in Eq. ^ . Inside the 
gaps between adjacent LLs of Eq. ([T]), the number of 
states is of the order of a. 



D. The Z2 nature of the topological properties 

Due to their connection to the 2D reduced version of 
the LL Hamiltonian, we still denote the low-energy bands 
of Eq. ^ as 2D parity breaking LLs. As shown in Eq. 

although these LLs are not exactly flat, their disper- 
sion over jz is strongly suppressed by the large value of 
a. If the chemical potential p, lies in the middle of the 
band gap, the Fermi angular momentum jkf,z is at the 
order of a. The classic radius of such a state is roughly 
It- As analyzed in Sec. Ill Bl for states with \ jz\ < a, two 
Hamiltonians Eqs. H]) and ([7]) share the same physics. 

Compared to the usual 2D LL states, the SO coupled 
LLs of Eq. ID) in the form of Eq. ^ are markedly 
different. The smallest length scale is not It, but the 
SO coupling length scale Iso = It 1 01 ^ It- Instead, we 
can use It as the cut-off of the sample size by imposing 
an open boundary condition at the radius of It- States 
with \ jz\ < a are considered as bulk states which localize 
within the region of p < It- States with \jz\ ~ a are edge 
states. 

We take the thermodynamic limit as follows. First, oj is 
fixed, which determines the LL gaps. Then we set m ^ 
and A — )■ 00 while keeping Igo = h/{mX) unchanged, such 

that It = \ — > 00. The number of bulk states scales 

linearly with a, and the level spacing scales as 1/a — 
at the Fermi angular momentum jkj,z ~ ct- 

The next important question is the stability of the gap- 
less edge modes. This situation is different from the usual 
2D LL problem, in which inside each LL for each value 
of angular momentum m, there is only one state. Those 
edge modes are chiral and, thus, robust against external 
perturbations. Since Eq. ^ is TR symmetric, for each 
filled LL there is always a pair of degenerate edge modes 
V'rir.iij on the Fermi energy, where rir is the LL index. 
Nevertheless, these two states are Kramer pairs under the 
TR transformation satisfying = — 1. In other words, 
the edge modes are helical rather than chiral. 

We generalize the reasoning in Ref. and[3 for topolog- 
ical insulators with good quantum numbers of lattice mo- 
menta to our case with angular momentum good quan- 
tum numbers. Any TR-invariant perturbation cannot 
mix these two states to open a gap. In other words, the 
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mixing term, 

Hmx = 5(V'|_D,„,.J^'02D,ri, -i, + h.c), (12) 

is forbidden by TR symmetry. On the other hand, if two 
LLs with indices and n'^ cut the Fermi energy, the 
mixing term, 

+ h.c), (13) 

is ahowed by TR symmetry and opens the gap. Conse- 
quently, the topological nature of such a system is char- 
acterized by the Z2 index, even though it is not clear how 
to define the Pfaffian-like formula for it due to the lack of 
translational symmetry.— Similarly to the 2D topological 
insulators based on lattice Bloch-wave states, in our case, 
if odd numbers of LLs are filled such that there are odd 
numbers of helical edge modes, the gapless edge modes 
are robust. 

Imagining an open boundary at p w It, we derive an 
effective edge Hamiltonian for these helical edge modes. 
As \jz\ a and taking the limit of a — >■ -l-oo, these 
edge modes are pushed to the boundary. We expand the 
spectra around jzjm- The edge Hamiltonian in the basis 
of jz can be written as 

Hedge = " V'". (W) 

jz 

where fj, = ^j^p-jzjm- The edge modes ipnr-,±jz around 
jzjm can also be expanded as 

V'n.j, = ^ "*(m+l)0 > V-n^.-i, = ri/;„^j^_.(15) 

\ yn^^ J 
fn^ and are real numbers parameterized as 

in,. = COS — , 5„,, = sm — , (16) 

which are determined by the details of the edge. We 
neglect their dependence on \jz\ for states close enough 
to the Fermi energy. The effective edge Hamiltonian can 
also be expressed in the plane-wave basis if we locally 
treat the edge as flat 

Hn^^edge = Vf (^Sm6„^[{p X 6^) • z)]{a ■ ir) 

+ cos ein^ (p X e^) • cr^^ - /i, (17) 

where is the local normal direction on the circular 
edge; both terms are allowed by rotational symmetry, 
TR symmetry, and the vertical mirror symmetry in such 
a system. Each edge channel is a branch of helical one- 
dimensional Dirac fermion modes. 

Equation ([T]) can be defined on the compact 5^ sphere, 
which takes the simple form 

H=^-u;L-a. (18) 



The eigenvalues of L • c? take Ih and — {l + l)h for the pos- 
itive and negative helicities of j± = / ± ^ , respectively. 
For convenience, we choose the parameter value of luj/h 
as a large half-integer, then for the lower energy branch, 
the energy minimum takes place at jo,+ = + 5 — Ii^/fi. 
The lowest LLs become SO-coupled harmonics with j+ — 
jo.+ and (2^0 -I- 2)-fold degeneracy. The gap between the 
lowest LLs and higher LLs is A = h^/{21), which is in- 
dependent of w. To take the thermodynamic limit, we 
keep / constant while increasing the sphere radius R, 
and maintain uj scaling with , such that the density of 
states on the sphere is a constant. 



III. THREE-DIMENSIONAL SPIN-ORBIT a p 
COUPLING IN THE HARMONIC TRAP 

In this section, we generalize the results in Sec. |ll] to 
three dimensions. We consider the ct -p-type SO coupling 
combined with a 3D harmonic potential 

HsD = -^^^ + ^mwV-A(-zfiV-c7). (19) 
2m 2 

Equation (IT9t possesses the 3D rotational symmetry, and 
TR symmetry of fermions with = —1. The parity 
symmetry is broken by the a ■ k term, and there is no 
mirror plane symmetry either. The quantities /go, It, ci, 
and fco are defined in the same way as in Sec. [Hi 

Although it is difficult to realize strong SO coupling in 
the form of u -p in solid-state systems, it can be designed 
through light-atom interactions in ultra cold atom sys- 
tems. We present an experimental scheme to realize Eq. 
(Ull) in Sec. |IVl 



A. Energy spectra 

Again, we consider the limit of strong SO coupling, i.e., 
a ^ 1. It is straightforward to generalize the momen- 
tum space picture in Sec. |IT]to the 3D case as presented 
in Ref. 46 and summarized below. The helicity opera- 
tor (T • fc is employed to define the helicity eigenstates of 
plane waves (t? • fc)'0j_|_ = itpj: Only positive- helicity 
states ipj:^ are kept in the low energy Hilbert space. The 
harmonic potential becomes the Laplacian operator in 
momentum space, and thus is equivalent to a quantum 
rotor subject to the Berry phase in momentum space as 
Vfp = im(zVfe — Ak)'^ . The moment of inertial is again 
/ — MkkQ and Mk = tf/{muj^). The Berry connec- 
tion Ak = iiipk ±l^fclV'fc ±) is the vector potential of the 
U(l) magnetic monopole. As a result, the angular mo- 
mentum quantization changes to that j takes half-integer 
values starting from i. The energy dispersion becomes 

EagiU) = fi^jU + l)/2/ = UU + l)/2a^)Etp, and each 
level is (2j -I- l)-fold degenerate. The radial quantiza- 
tion is the same as before. Thus the dispersion can be 
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summarized as 



2 2a2 



(20) 



where rir is the band index, or, the LL index, and j is 
the angular momentum quantum number. 



Dimensional reduction from the 4D Landau 
level Hamiltonian 



Following the same logic as in Sec. IIIBl we present the 
dimensional reduction from the 4D LL Hamiltonian [Eq. 
([22]) ] to arrive at a 3D SO coupled Hamiltonian closely re- 
lated to Eq. (HH). The 3D LL Hamiltonian, Eq. ©, can 
be easily generalized to arbitrary dimensions by combin- 
ing the n-D harmonic potential and the n-D SO coupling 
between orbital angular momenta and fermion spins in 
the fundamental spinor representations.— In four dimen- 
sons, there are two non-equivalent fundamental spinors, 
both of which have two components. Without loss of 
generality, we choose one of them as 



(21) 



where i,i = 1,2, and 3. The orbital angular momentum 
operators are defined as Lij = —ihxi\/j + ihxjS/i where 
i,j = 1,2,3, and 4. The 4D LL Hamiltonian in the flat 
space is deflned as 



AD.LL 



2m 



4 

E 

- E 

l<i<j"<4 



E^^ 



(22) 



which possesses TR and parity symmetry. 

The ^th-order 4D orbital spherical harmonics coupled 
to the fundamental spinor can be decomposed into the 
4D SO-coupled spherical harmonics in the positive- and 
negative-helicity sectors, where Lijatj take eigenvalues of 
Ih and —{I -f 2)ft, respectively. The eigen wave functions 
of Eq. (|22l) in the positive-helicity channel are disper- 
sionless with respect to / as En^,+ = {2nr + 2)hu!. Their 
radial wave functions are i?„^_/(r) = r'e^'' ^^^t p(^—n^^ 
2,r^//|,), where F is the standard confluent hypergeo- 
metric function. With an open boundary of an 5*^ sphere, 
each filled LL contributes to a gapless surface mode of 3D 
Weyl fermions as 



H: 



3D .surf ace — '^f^r,i^ijPj f-^i 



(23) 



where is the unit vector normal to the sphere. The 
topological index for such a 4D LL systems with TR sym- 
metry is Z rather than Z2. 

We perform the dimensional reduction on Eq. ([^^ 
from four to three dimensions. We cut a 3D off-center 
hyper-plane perpendicular to the fourth axis with the 



interception X4 = wq Within this 3D hyper-plane of 
(xi, X2, 2:3, X4 = Wo), Eq. reduces to 



H 



3D, rede 



= H- 



3D.SO 



— LUL ■ (T, 



(24) 



where the first term is just Eq. (|22p with the SO- 
coupling strength A — lowq. It contains another SO- 
coupling term, L ■ a, and its coefficient is the same as the 
harmonic trapping frequency. Similarly to the previous 
reduction from three to two dimensions, here we have 
a — It/Iso — \wo\/It- At wq — 0, Eq. (IMl) becomes the 
3D LL Hamiltonian of Eq. ([5]) with parity symmetry. If 
Wq 7^ 0, the (T-p term breaks parity symmetry. Following 
the same reasoning as in Sec. Ill Bl Eqs. ([T^ and ([M)) 
share the same physics for eigenstates with j < a in the 
case of a ^ 1 . 

Similarly as before, we construct an off-center solution 
to the 4D LL problem. We use r to denote a point in the 
subspace of a:i_2,3, and f2 as an arbitrary unit vector in the 
X1-X2-X3 space. We consider the plane of il-x^ spanned 
by the orthogonal vectors Cl and X4. It is easy to check 
that the following wave functions, which depends only on 
coordinates in the Q-X4 plane are the lowest LL solutions 
to the 4D LL Hamiltonian, Eq. 



{f-n + ix4)e ^'t a, 
where a^^ = (cos |, sin |e"'^)'^ satisfies 



(25) 



(26) 



In this set of wavef unctions, both the orbital angular mo- 
mentum and spin are conserved and added up; they are 
called the highest weight states in group theory. In fact, 
these states can be rotated into any plane accompanied 
by a simultaneous rotation in the spin channel. Based on 
the structure of the highest weight states, we can still de- 
fine the magnetic translation operator in the Cl-Xi plane 
along the X4 axis as 



'^0x4(^0^4) =exp 



(27) 



Applying this operator to the Gaussian pocket of the 
solution with / = in Eq. ([25)) . we arrive at the off- 
center solution 



9;^ 



(28) 



This solution, however, breaks the rotational symmetry. 
In order to restore the 3D rotational symmetry around 
the new center (0,0,0,100), we perform a Fourier trans- 
formation over the direction of VL as 

(29) 

M -I- ^ . Please note that due to 
the singularity of an over the direction of fi, monopole 



where j = I + i and jz 
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spherical harmonics, i 1 ,„+ 1 (rJ), are used instead 
of regular spherical harmonics. 

Again, noting that dx^ipio-jj, {r, X4)\xi,wo = 0, we sim- 
ply set Xi — Wo', then it is simple to check that the re- 
duced 3D wave functions 

^3D j j, (r) = i^4D-j,j, {r, Wq) (30) 
are the solutions to Eq. (l24l) for the lowest LLs as 

(31) 

ipsD.jj^if) can be simplified as 

J, 2 

X r_,,-,+i,,,(f^r)}, (32) 

where fco = wo/l'^ = mX/h and A = wquj; ji is the Zth- 
order spherical Bessel function. Y±jjj^^s are the SO- 
coupled spherical harmonics defined as 

with a positive eigenvalue of ^/i for a ■ L, and 

with a negative eigenvalue of — (/ + l)h for a ■ L. 

The difference between Eq. and Eq. ([19]) is the 

term cf • L, whose effect is weakened as the distance from 
center r gets small. The radial distributions of ji{kor) in 
Eq. and Jm(/cop) in Eq- ® are similar. Following 
the same reasoning presented in Sec. Ill Bl in the limit of 
a 1, we can divide the lowest LL states of Eq. ([5^ 
into three regimes: j < a, j ^ , and a < j < a^. At 
j < a, the classic orbit radius scales as rc.j ~ ^It < It- 
Again in this regime, the effect of (? • L is a perturbation 
of the order of rcj^/zo ^ 1; thus the two Hamiltonians, 
Eq. ([M]) and Eq. p^ . share the same physics. Similarly, 
in the regime of j ^ a^, a ■ L dominates, and the physics 
of Eq. ([24]) conies back to the 3D LL Hamiltonian, Eq. 
dS]), while that of Eq. ([TO)) is no longer LL-like. 

C. The Z2 helical surface states 

Following the same reasoning as in Sec. IIID| we denote 
the low-energy bands of Eq. (fT9|) as 3D parity breaking 
LLs. For the lowest LL, below the energy of the bottom 
of the second LL, the angular momentum j takes values 
from i to the order of a at which the radius of the LL 
approaches It- For this regime j < a, Eqs. and ([M)) 
share the same physics. Again, the smallest length scale 



is the SO coupling length scale Iso = Ir/ct ^ It- States 
with \j\ ^ a are considered bulk states which localize 
within the region p ^ It- States with \jz\ ^ a are edge 
states. The number of bulk states scales linearly with c? . 

Now we impose an open boundary condition of an S*^ 
sphere with radius r k, It, and consider the stability of 
the edge modes against TR invariant perturbations. Let 
us consider one filled LL. The Fermi energy lies between 
the gap, and thus cuts the dispersion at surface states. 
In the limit of a 00, the energy level spacing between 
adjacent angular momenta j and j -t- 1 scales as huj /a — > 
for surface modes with j ^ a. Thus we can always 
choose the Fermi angular momentum jf satisfying jf = 
21 + ^. For this value of jf, there is an odd number of 
2Z -I- 1 Kramer pairs between ipjf.±jz for jz = i to jf. 
Again according to the reasoning of the Z2-classification 
in Refs. 6 and 7], these states cannot be fully gapped out 
by applying TR invariant perturbations. Certainly, for 
those states with j = 21 + ^ close to the Fermi energy, 
they can be fully gapped, but they are only part of the 
spectra, and do not change the topological properties. 
Again, if two LLs with different indices rir and n'^. cut the 
Fermi energy, the zero energy states at the Fermi level 
can be fully gapped out. Thus, the topological nature of 
Eq. (Ull) is Z2. 

We further present the effective surface Hamiltonian 
for surface modes in the limit of j/ ~ a — -l-oo. The 
effective surface Hamiltonian of the 3D topological in- 
sulators with the spherical boundary condition has also 
been discussed in Refs. [sil andlH. The surface Hamilto- 
nian in the eigen-basis of j and jz can be written as 

= E( v'^'' " (33) 

where fi = -j^j/. The construction of the accurate 
surface Hamiltonian in the plane-wave basis depends on 
the detailed information of surface modes ipjjz (^1 ^r) for 
j sa jf and, thus, is cumbersome. Nevertheless based on 
the symmetry analysis, we can write the general form as 

H„^,edge = Vf^^Sm9n,.iP X (?) • 6^ 

-I- cos6'„^ [p- a — ip- er){(T ■ e^)] | — A^.(34) 

where is the local norm direction on the S'^-sphere. 
Both terms obey the local SO{2) rotational symmetry 
around the and TR symmetry. The first Rashba term 
also obeys the vertical mirror symmetry, while the second 
term does not. The second term favors the spin aligning 
with the momentum, while the second favors a relative 
angle of 90°. For a general value of 9n^, which is deter- 
mined by the non-universal surface properties and 0„^, , 
Eq. (I34p determines a relative rotation between spin and 
momentum orientation at the angle of It is still a 
helical Dirac Fermi surface. 
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4^ 



FIG. 2: (Color online) Level diagram for atom-laser coupling. 
Four lower energy levels are coupled to two excited levels to 
compose a hybrid tripod and tetrapod configuration. 



B,\+ b) 



42 



IV. EXPERIMENTAL REALIZATION FOR 3D 
SO COUPLING 



FIG. 3: Energy levels of the atom-laser coupling Hamiltonian 

Eq. uni 



In the ultra cold atom context, there has been great 
progress in the synthetic gauge field, or, artificial SO cou- 
pling from light-atom interactions^"^. Experimentally, ar- 
tificial SO coupling has been generate in ultra cold atom 
systems 4^ Two dimensional Rashba and Dresselhaus SO 
coupling in the harmonic potential has been proposed 
using a double-tripod configuration*^^ Since the pseudo- 
spin degrees of freedom are represented by the two lowest 
energy levels, this scheme is immune to decay due to col- 
lision and spontaneous emission process.— 

In this section, we propose the experimental realization 
for the 3D SO coupling of the g-p type in Eq. ([TO]). Here 
we generalize the scheme in Ref. |5J to a combined tripod 
and tetrapod level configuration as depicted in Fig. [2] 
Three internal levels |1), |2), and |3) couple the excited 
state I a) to form a tripod configuration. A tetrapod- 
like coupling is formed by coupling the four levels |1) — 
1 4) to the common excited state \h). The single-particle 
Hamiltonian reads 

H = ^ + ]^muj\-' + Hal, (35) 

where m is the mass of the atom; Hai represents the 
atom-laser coupling. In the interaction picture, Hai can 
be written under the rotating wave approximation as 



To further simplify the model, here we assume {Ba\Bi,) = 
0, which can be achieved by choosing 



(j = l,2,3;m = a,6)(39) 



with 01 



.ja - [i - 2)27r/3 and Ojb = -(j - 2)27r/3. We 
also choose 0.^, — ri4e*^'^'' '^+^*\ and set Q.c = Via, = 
51c cos 0, and = sin (j). Using these notations, Hai 
is simplified as 



Ha 



n,{\a){Ba\ + \b){Bb\) + h.c. 



(40) 



where \Bb) = cos (j)\Bb) + s'm(f)\l) and |4) = 
g-i(/c4 r+e4)|^^^ The above Hamiltonian supports three 
pairs of degenerated eigenstates with energy difference 
h\Qc\, as depicted in Fig. Explicitly, the eigen- 

vectors are written as 

\Ba) + \a) \Bb) + \b) 



\Gi) = ^^^1^: IG2) = 
IG3) - \D), 
IG3) = \Ge 



V2 



IG4) = \B^), 
\Bb) 



\b) 



V2 



V2 



(41) 



Hal = E {^iM{l\ + ^2rn\m){2\ + n„^\m){3\ whereJi?)=E,e-'=^-1j)/V3and|S,^)=sin^|B,)- 



h.c.} -n[n4b\b){4\ + h., 



(36) 



where flim are the corresponding Rabi frequencies be- 
tween the internal states \i) and |m) with m — a,b . 
We introduce the following two bright states 



\B„,) = {nUi) + nu 



(37) 



cos 0|4) . Therefore, the two degenerate ground states can 
be used as pseudo-spin 1/2 degrees of freedom. 

If the trapping frequency satisfies w ^ \^c\, according 
to the adiabatic approximation, we neglect the coupling 
between the ground-state manifold and other states. 
Therefore, atoms in the subspace spanned by |Gi) and 
IG2) evolve according to the effective Hamiltonian 



where m = a, 6 and Qm = + l^2mp + jOsmp. 

The atom-laser coupling can be rewritten as 

Hal = ~h{na\a){Ba\+h.C.} 

- hi^nb\b){Bb\ + nib\b){M + h.c.y (38) 



H, 



jp- 
2m 



+ -muj r + <i>, 



(42) 



where the non-Abelian gauge potential A is a 2 x 2 matrix 
with the elements 



A,,^ih{G,\V\Gj 



(43) 
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F = U 




(Color online) Energy level scheme for alkali atoms 

9 



The Zeeman sublevels of two hyperfine states F — ^ 



FIG. 

and F = ^ can be used to fulfill our requirements. Lines 
or curves with an arrow indicate effective transitions between 
different magnetic levels which can be implemented using res- 
onant Raman processes. Other levels, which are not involved 
in the scheme, are not shown. 



where i,j — 1,2; $ is a scalar potential induced by the 
coupling laser beams. 

An isotropic 3D a • p-like SO coupling can be obtained 
by a 3D set-up of laser configurations as 



fci = k{- 



1 



V3 



0), = Ai(0,l,0), 



In this case, the corresponding vector and scale potential 
are calculated as 



— = 0.166K[aj:e^j: + ayCy + {a^ - I)e^] , 



0.445^/. 

2m 



(45) 



The $ term is a constant and, thus, can be dropped off. 
The Abelian part in the gauge potential Az is a con- 
stant, which can be absorbed by a gauge transformation. 
Consequently, the remaining constant non- Abelian gauge 
potential behaves as a tr • p type SO coupling. 

The above-considered level structure can be found for 
example, in alkali atoms with large spins. Figure S] shows 
the hyperfine ground state manifolds of for '^°K 

atoms under an external magnetic field. The energy lev- 
els |1) ^ |4), |a), and \b) can be selected as different Zee- 
man sublevels of = | and -F" = | ■ Using the notation 
of \FMp) to denote each state, we choose |1) = ||,— 
|2) = |3) = |4) = \a) = 



and 1 6) = The coupling between different levels 

is achieved for example, by using two laser beams under 
second-order resonant Raman process. The two lasers 
can be chosen to be circularly polarized and tt polarized, 
respectively, in order to satisfy the selection rule. Fi- 
nally, wave vectors of individual laser beams can also be 
adjusted so that Eq. (f45|) is fulfilled. 



V. CONCLUSION AND OUTLOOK 

We have studied rotationally and TR symmetric LL 
systems in both 2D and 3D systems with breaking par- 
ity symmetry, whose topological properties are character- 
ized by the Z2 class. These Hamiltonians are simply 2D 
harmonic potentials plus Rashba SO coupling, or 3D har- 
monic potentials plus tr-p-type SO coupling with a strong 
SO coupling strength. For low-energy bands, the disper- 
sions over angular momenta are strongly suppressed by 
SO coupling, to be nearly flat. Up to a small difference 
which can be treated perturbatively, these Hamiltonians 
can be systematically investigated through dimensional 
reduction on the high-dimensional LL problems by cut- 
ting an off-center plane in the 3D LL Hamiltonian or an 
off-center hyper plane in the 4D LL Hamiltonian. The 
parity breaking LL wavefunctions in two and three di- 
mensions are presented explicitly. With open boundary 
conditions, helical edge states are found in two dimen- 
sions, and surface states are found in three dimensions. 
These states can be realized in ultra cold atom systems 
in a harmonic trap combined with synthetic gauge fields, 
i.e., artificial SO coupling. In particular, we propose an 
experimental scheme to realize the 3D Hamiltonian. 

The above dimensional procedure can be straight- 
forwardly generalized to arbitrary dimensions based 
on our previous construction of high dimensional LL 
Hamiltonians,— and so can the general parity breaking 
LL wavefunctions in N dimensions. The nice analytical 
properties of the 2D and 3D LL wave functions breaking 
parity symmetry also provide a good opportunity to fur- 
ther construct many-body wave functions of the factional 
topological states. These properties will be investigated 
in a future publication. 
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